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Repulsively bound pairs of particles in a lattice governed by the Bose-Hubbard model can form stable incom- 
pressible clusters of dimers corresponding to finite-size n — 2 Mott insulators. Here we study the dynamics of 
hole defects in such clusters corresponding to unpaired particles which can resonantly tunnel out of the cluster 
into the lattice vacuum. Due to bosonic statistics, the unpaired particles have different effective mass inside 
and outside the cluster, and "evaporation" of hole defects from the cluster boundaries is possible only when 
their quasi-momenta are within a certain transmission range. We show that quasi-thermalization of hole de- 
fects occurs in the presence of catalyzing particle defects which thereby purify the Mott insulating clusters. We 
study the dynamics of one-dimensional system using analytical techniques and numerically exact t-DMRG sim- 
ulations. We derive an effective strong-interaction model that enables simulations of the system dynamics for 
much longer times. We also discuss a more general case of two bosonic species which reduces to the fermionic 
Hubbard model in the strong interaction limit. 

PACS numbers: 37.10.Jk, 03.75.Lm, 67.80.dj, 05.30.Jp, 



I. INTRODUCTION 

Quantum particles in lattice potentials, e.g. electrons in 
crystals, have been studied since the early days of quantum 
theory |1, 2|. With the development of artificial (optical) 
lattice potentials for cold neutral atoms |3|, bosonic lattice 
models are recently attracting increased interest as well (H, 
with the Bose-Hubbard model (BHM) |5 1 being an important 
example. A remarkable phenomenon entailed by the BHM 
is that pairs of strongly interacting bosons can form tightly 
bound "dimers" both for attractive and repulsive interactions 
|6~9|. In free space, or in the presence of energy dissipation, 
the repulsive interaction inevitably leads to pair dissociation. 
In a lattice, however, the kinetic energy of each particle is 
restricted to the values in the allowed Bloch band. Conse- 
quently, two co-localized particles in a dissipation-free lattice 
remain tightly bound together as a dimer when their interac- 
tion energy U exceeds the kinetic energy of free particles ~ J 
within the Bloch band. 

In a previous publication |9 |, we have studied the many- 
body dynamics of the repulsively-bound dimers of bosons. 
Due to virtual transitions of the dimer constituent particles, 
the dimers at the neighboring lattice sites strongly attract 
each other, with the corresponding interaction energy exceed- 
ing the dimer tunneling energy by a factor of 4. For many 
dimers on the lattice, it is then energetically favorable to 
form dynamically stable "droplets", constituting incompress- 
ible Mott-insulating (MI) clusters with the number of parti- 
cles per site of exactly n = 2. Inevitable imperfections in the 
preparation process would typically cause such MI clusters to 
contain hole and particle defects corresponding, respectively, 
to unpaired and excess particles (monomers and trimers). An 
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important question is thus how to purify the system of the de- 
fects reducing thereby the entropy. In the present paper, we 
discuss a mechanism of self-purification of stable MI clus- 
ters of dimers surrounded by lattice vacuum. We study the 
dynamics of defects in one-dimensional system by analytical 
calculations and numerical many-body simulations. 

Within the cluster, hole and particle defects can propagate 
via resonant single-particle hopping with enhanced amplitude, 
which stems from the bosonic statistics of the surrounding 
n = 2 MI environment. Outside the cluster, hole defects 
correspond to free particles. Since their tunneling energy J 
is much larger than the monomer-dimer interaction energy 
~ J'^ /U fTol, hole defects are not bound to the cluster and 
can "evaporate". However, the widths of the single particle 
Bloch band is twice larger inside the cluster than outside of 
it, therefore only the hole defects with energies in the center 
of the band can penetrate the cluster boundaries and evapo- 
rate into the lattice vacuum, while in the absence of quasi- 
momentum redistribution, low- and high-energy hole defects 
will remain in the cluster. We show that the presence of par- 
ticle defects leads to efficient "thermalization" of the hole de- 
fects via quasi-momentum redistributing collisions. Hence, 
very few such "catalyzing" particle defects can purify the MI 
cluster. 

Before continuing, we note a recent relevant work fTTIl 
dealing with fermionic dimers described by the Hubbard 
model. After preparing a cold atomic gas with filling of n 2 
in the trap center, followed by turning off the trap, the hole de- 
fects will simply tunnel out of the cluster into the vacuum. For 
fermions, however, the remaining cluster is not stable, since 
the effective second-order tunneling of the on-site pairs is not 
restrained by the interaction between the pairs. 

Figure [T] illustrates the main physics studied in this paper, 
which is organized as follows. In Sec.[ll|we outline the prop- 
erties of the pure dimer clusters 0. We then introduce in 
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FIG. 1 : (Color online) Physical models studies in this paper, (a) The 
Bose-Hubbard model, Eq. ([TJ. (b) Monomer (hole defect) effec- 
tive hopping, Sec. |lll| (c) Single defect effective theory. Sec. |lll| 
(d) Trimer (particle defect) effective hopping, Se c.[rv| (e) Single hole 
defect in two species Bose-Hubbard model, Sec.|VI| 



(hole defect) from a domain wall separating the dimer cluster 
and the vacuum. The quasi-momentum redistribution of a hole 
defect upon collisions with a particle defect in the lattice with 
periodic and open boundaries is studied in Sec.|lV] In Sec.[V| 
we present the results of many -body numerical simulations for 
a realistic system with several hole and particle defects in a 
dimer cluster surrounded by empty lattice. Finally, in Sec. VI 
we discuss the case of two bosonic species, which is more 
flexible theoretically, but is demanding experimentally. In the 
limit of infinite intra-species interaction, it contains the special 
case of the Hubbard model |11|, since in one dimension and 
in the absence of double occupancy, bosons and fermions are 
equivalent through the Jordan- Wigner transformation. Much 
of the involved technical details are deferred to Appendices [A| 
|Bl|C]and|Dl 



II. REPULSIVELY BOUND DIMERS 

The underlying Hamiltonian for our system is that of the 
BHMlll 

H = -J^(6]6,+i +H.a.) + '^U^bybybfb,, (1) 



where b and b^ are bosonic annihilation and creation opera- 
tors, J is the particle hopping rate between adjacent lattice 
sites j, j + 1, and U is the contact interaction between the 



particles on the same lattice site. Throughout this paper, we 
assume that the on-site interaction is the dominant energy pa- 
rameter, U ^ J. 

Considering first a lattice containing only zero or two par- 
ticles per site, we do not allow the dimer occupation number 
in the lattice to exceed unity. Adiabatically eliminating all the 
states with odd number of particles per site, we obtain for the 
dimer s an effective Hamiltonian |9 | that contains only terms 
with characteristic energies on the scale of J'^/U <C J: 

H = -J^(4c,+i +H.a.) + ^^44+iC,+ic„ (2) 



where J = — 2J^/[/ is the dimer hopping rate and B = 
— 16J^/U is the nearest neighbor interaction. The dimer cre- 
ation Cj and annihilation cj operators satisfy the hard-core bo- 
son commutation relations 



i=j: {cj.Cj} = 0, {cj,c]} = 1. 



(3a) 
(3b) 



Hamiltonian ^ can be mapped onto that for the spin-| XX Z 
model (UllTSl with the anisotropy parameter A = B/2J = 
4. For A > 1, we are in the ferromagnetic, Ising-like regime, 
and a cluster of dimers, corresponding to a lattice domain with 
maximum magnetization, is dynamically stable. To under- 
stand this in terms of dimers, observe that, for any U{:^ J), 
the maximal kinetic energy 2 J gained by releasing a dimer 
from the cluster boundary is small compared to the binding 
energy B of the dimer to the cluster. 

The stability of the dimer cluster is an intrinsic feature of 
the BHM. It is rooted in the bosonic amplification of the inter- 
site hopping of the particles, which in turn enhances the ef- 
fective (second-order) nearest neighbor interaction B. For the 
fermionic Hubbard model discussed in |11| in the context of 
defect evaporation from a dimer cluster, i.e., a band insulator, 
that A = 1 (5 = 2 J), which means that 



VI 



we show in Sec. 
the cluster is unstable and the dimers will diffuse away. 



III. SINGLE DEFECT MODEL IN THE 
STRONG-INTERACTION LIMIT 

The dynamics of dimers is rather slow, as it is governed 
by the small characteristic energies ~ J'^/U, but the dy- 
namics of monomers is much faster, involving large single- 
particle hopping rate J. We can thus retain only the con- 
tributions on the scale of J, which results in a very sim- 
ple and transparent effective theory for the monomers. For 
a monomer in the cluster (hole defect), the bosonic statistics 
plays an important role: it increases the hopping amplitude of 
the monomer in the environment of dimers by a factor of 2, see 
Fig.[TJb). As a result, the kinetic energy of the monomer in the 
dimer cluster is Ej^ = — 4Jcos(/c), while in the vacuum it is 
Ek = — 2Jcos(/c), where k G [— tt, tt] is the monomer quasi- 
momentum quantified by the phase change between neighbor- 
ing lattice sites. Therefore the monomer will be confined to 
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FIG. 2: (Color online) Transmission probability T{k) [cf. Eq. ( |A9| l] 
for various a = Jb/ Ja- 

the cluster if its quasi-momentum is not inside the transmis- 
sion region given by 

k e (-27r/3, -7r/3) U (tt/S, 27r/3), (4) 

up to a correction due to small interactions of the order of 
J'^/U <C J which we neglected. 

Consider the scattering of a monomer from the domain wall 
between the dimer cluster occupying sites j < (region A) 
and the vacuum at sites j > (region B), see Fig.[TJc). The 
local bare particle number is rij = 2 for j < and rij = for 
j > 0. The particle number at j = depends on the position 
i of the monomer = 1: inside the cluster z < we have 
no = 2, at the boundary i = obviously no = 1, while out- 
side the cluster (in the vacuum) i > leads to no = 0. Hence 
the position of the wall shifts upon the monomer crossing the 
boundary, which should be taken into account when consider- 
ing many defects. The hopping rate of the monomer is J a for 
sites j < and Jb for sites j > 0. The effective Hamiltonian 
for a single monomer then reads 

H = - Ja ^(a]a^-+i +H.a.) - Jb +H.a.), (5) 

j<0 j>0 

with Ja = 2 J and Jb = J. 

In Appendix |A| we calculate the exact transmission proba- 
bility T{k) of a particle crossing a domain wall in a system 
described by Hamiltonian ([5]). The results are illustrated in 
Fig. [2] for various a = Jb/ Ja- The values of a = 1/2 and 
a = 2 correspond, respectively, to the single particle leaving 
the dimer cluster and entering it from vacuum. 

IV. QUASI-MOMENTUM REDISTRIBUTION OF THE 
DEFECTS 

We have seen above that a hole defect can leave the MI 
cluster only if its quasi-momentum is within the transmission 
range of Eq. (|4]), while a defect with the quasi-momentum 



outside the transmission range will remain in the cluster in- 
definitely. Hence, to completely purify the cluster of hole 
defects, their quasi-momenta should be continuously redis- 
tributed over the entire range of k G [— tt, tt]. In two or more 
dimensions, collisions between identical particles can redis- 
tribute the absolute values of their quasi-momenta, and we 
therefore expect the evaporation of all the defects through the 
cluster boundaries after a few collisions. In a one dimensional 
lattice, however, collisions of two particles interacting via any 
finite range potential can only exchange their quasi-momenta 
or leave the quasi-momenta unchanged 17] |8l [141 El (for in- 
distinguishable particles both outcomes are equivalent), pro- 
vided that the lattice is deep enough so that the large band gap 
precludes interband transitions. Similarly, collisions with the 
fixed boundaries can only reverse the quasi-momentum of a 
particle. The simplest quasi-momentum redistribution mech- 
anism is then three particle collisions. This happens at a rate 
proportional to the defect density squared, which is too slow 
for cold atom experiments. 

In the dimer cluster, in addition to the hole defects 
(monomers), we may have particle defects (trimers) with dif- 
ferent effective mass. The hopping rates of a monomer and 
a trimer in the cluster are J^ = 2 J and Jt = 3 J, respec- 
tively. Fig. [ijd). Before collision, their quasi-momenta are ka 
and kt, while conservations of quasi-momentum, ka -\- h = 
k'^ + k[, and energy, Ja cos{ka) + Jt cos(/ct) = Ja cos(/c^) + 
JtCos(/c^), during the collision determine the new quasi- 
momenta /c^ and k[ via 

Ja COs{ka) + Jt COs{kt) = Ja COs(/c^) + Jt COs{ka h - k'^) . 

(6) 

If there is a collision with the wall, or a third defect of either 
kind, before this process is reversed, all energetically allowed 
combinations of ka , kt can be assumed, as will be verified be- 
low by exact numerical simulations. 

A. Two classical particles 

The timescale for quasi-momentum redistribution can be 
calculated from purely classical considerations. A monomer 
or a trimer moving in the MI cluster has a kinetic energy of 
Ek^ = — 2 J^ cos{k^) and the corresponding group velocity 
of = 2 J^ sin(A:^) [/i = a, t]. 

Consider first two wave packets in a periodic lattice of 
length L. After a collision (the defects can not penetrate each 
other), their velocities are assumed to be Va = 2Ja sin{ka) < 
Vt = 2Jt sin(A:^). The next collision happens after time 



2 2 Jt sm{kt) - Ja sin(A:a) ' 

and the new quasi-momenta are determined by Eq. ([6]). It fol- 
lows that the time interval between all subsequent collisions 
is the same tc/2, since Eq. ^ and 

L-1 1 

2 Jtsm{kt) - Jasm{ka) 
L-1 1 

2 -Jt Sm{ka ^h- k'a) + Ja ^^"^{k'a) 
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FIG. 3: (Color online) Dynamics of the quasi-momentum distribu- 
tion for the monomer (left column) and trimer (right column) in a 
lattice of L = 64 sites. The initial quasi-momenta are ka = ^tv and 
kt — — ^TT. Upper panels correspond to periodic boundary condi- 
tions, where the markers on the right indicate multiples of the revival 
time tc ~ 46.63/Jt. Lower Panels are obtained for open boundary 
conditions. Dashed vertical lines mark the transmission regions for 
monomer quasi-momenta as per Eq. (|4]). 

always have a common solution. 

In the presence of a wall, or a third defect, the quasi- 
momenta can take any values energetically allowed. A re- 
vival is not expected, but now is an effective rate of quasi- 
momentum redistribution. It is essentially given by J over the 
mean free path, i.e., it is proportional to J times the average 
defect density, which is indeed much faster than the rate of 
three particle collisions. 



B. Two quantum particles: Numerical simulations 

We simulate the quantum dynamics of the hole and particle 
defects in a dimer cluster using the two-particle Hamiltonian 
in quasi-momentum space, see Appendix |B] Each defect is 
initially prepared in a quasi-momentum eigenstate, with the 
combined state given by 

1 ^ 

where creates a monomer and a trimer at site j of a finite 
lattice filled with dimers, playing here the role of an effective 
vacuum |vac). 

Eigure [3] shows the results of our numerical simulations. 
In the case of periodic boundary conditions, the dynamics is 
mainly classical; only two values of quasi-momentum k are 
assumed by each particle, and after the classical revival time 
tc, the quasi-momentum distribution is restored to the initial. 
Eor open boundary conditions, however, we observe fast re- 
distribution of quasi-momenta, and already the first revival 
is hardly noticeable. We may therefore conclude that a sin- 
gle trimer can catalyze the redistribution of quasi-momenta of 



monomers, making the evaporation of almost all hole defects 
possible, provided that their average kinetic energy is initially 
close to the center of the band. This will be verified by the 
following many-body calculations. 



V. MANY-BODY NUMERICAL SIMULATIONS 

To study the dynamics of several defects under experimen- 
tally realistic conditions, we use a sufficiently long lattice 
that can accommodate dimer clusters spanning a few dozen 
sites. The complete Hilbert space for such a system is too 
large to be amenable to exact diagonalization treatments. We 
therefore resort to time-dependent density matrix renormal- 
ization group (t-DMRG) methods |[T6|[T7l, specifically, the 
time evolving block decimation (TEBD) algorithm 1 18] using 
the matrix product state (MPS) formalism. Even then, how- 
ever, simulating the full BHM is a difficult task. This is due 
to sizable quantum fluctuations present even in the pure dimer 
cluster for any finite interaction strengths U/ J. These fluctua- 
tions contribute to the many-body entanglement and consume 
much of the computational resources required to simulate the 
dynamics of the defects. We therefore introduce an effective 
model for the defects only. 



A. Many defect effective theory in the strong-interaction limit 

Since the states with different number of particles per site 
have energies separated by multiples oilKC^ J), the numbers 
of monomers, dimers, and trimers in a lattice are, to a good 
approximation, conserved separately. This allows us to treat 
the monomers, dimers, and trimers as distinguishable species, 
each represented by hard-core bosons, Eqs. ([3]). Eurthermore, 
as discussed in Sec. |Il| dimers forming stable clusters do not 
contribute to the dynamics of the system. Eor our initial con- 
ditions, typically containing a single cluster, we can thus re- 
formulate the problem as one of the hole and particle defects 
moving on the background of dimers or vacuum, with the spa- 
tial configuration of the dimer cluster entering the effective 
Hamiltonian for the defects only as a parameter. 

We define the reference system in which the pure dimer 
cluster occupies certain lattice sites while all the defects are 
placed at the beginning (left side) of the lattice. As the de- 
fects move in the lattice, the effective hopping rates depend 
on whether they are inside or outside the MI cluster. In turn, 
the position of the cluster depends on the positions of the de- 
fects, since each defect crossing the system from the left to 
the right shifts the position of the dimers, and the cluster as a 
whole, by one site to the left. The effective Hamiltonian for 
the defects can then be cast as 

j = l nr=0 j = l 

where P^\2 l] projector onto the subspace containing 

exactly hole and particle defects on sites j + 2 to L, while 
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each local operator i^j®^ acts on sites j and j 



1 as 



[O] 



jf + H.a.)%a,+ia]+ia]. (10) 



Here and (fj and ij) are the hard-core bosonic creation 
and annihilation operators for the monomers (trimers). The 
function 6(j) is initialized for all j with respect to the refer- 
ence system, and it can take two values: 6(j) = 1 for site 
j + 1 being empty (vacuum) and 6(j) = 2 for site j + 1 
containing a dimer. Then the hopping rates for the monomers 
are ji^^ = J and J^^ = 2 J, and for the trimers are = 
(they can not move on an empty lattice in first order in J) and 
Jfl = 3 J. 

Note that since the effective Hamiltonian ([TO]) contains 
two species of particles with hardcore interactions, it can not 
be mapped onto a model of free fermions via the Jordan- 
Wigner transformation (which is possible for identical hard- 
core bosons). The dynamics is therefore non-trivial and ac- 
tual calculations again require numerical many-body (TEBD) 
techniques. The practical advantage of the effective model — 
besides the largely reduced number of particles — is that the 
fast timescale /7~^ is eliminated from the system's dynamics 
and in our numerical simulations we can choose Trotter steps 
on the time scale ^ Further discussion on the effective 
defect model is given in Appendix [C] 



B. Initial states 

In our numerical calculations, we use several typical con- 
figurations of the defects in the lattice, each configuration de- 
scribed by a pure quantum state. Various coherent and inco- 
herent superpositions of such configurations would represent 
mixed initial states. 

We consider piecewise product states. A MI segment of 
length / contains fixed number of particles n at every site (n = 
2 inside the dimer cluster and n = in the vacuum). 



vac), 



(11) 



with I vac) denoting the true vacuum. Each segment can con- 
tain an additional defect. For a defect localized as site j, we 
use the notation 



I \n 



(n>l), 



(12a) 
(12b) 
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FIG. 4: (Color online) Density of monomers (left column) and 
trimers (right column) in the n — 2 MI cluster of 24 sites surrounded 
by empty lattice, |-)§2» on both sides. In the top panels, the ini- 
tial state of the cluster | — 7r/2)g~ |-)||7r/2)g~ corresponds to two 
monomers at the center of the band moving to the left and right. 



In the central panel, the initial state |7r) 



^|0)8 corresponds to 



two monomers at the upper and lower band edges. In the bottom 
panel, the initial state |7r)g~ |7r/2)g^|0)g~ is the same as in the cen- 
tral panels plus a particle defect at the center of the band, moving 
to the right. The interaction strength is [/ = 100 J. The density of 
monomers (trimers) corresponds to the probability of finding exactly 
one (three) particles at a given site. A TEBD 1 18 | algorithm with 
bond dimension x = 200 is used for the time evolution with a fourth 
order Trotter decomposition and time step size 1/50J, with particle 
number conservation explicitly included in the MPS 1.19 J . 



k, which must be a multiple of 27r//, as 



(13a) 



We prepare the cluster by joining MI segments with and 
without defects. Since we are only interested in low defect 
densities, we do not construct segments containing multiple 
defects. In order to perform TEBD simulations, the initial 
states have to be represented in the MPS form, which is dis- 
cussed in Appendix [D] 



C. Numerical results 



with ± corresponding, respectively, to a particle and a hole 
defect. Similarly, we denote a defect with quasi-momentum 



Figures [4] and [5] show the time evolution of defects in a 
n = 2 MI cluster surrounded by vacuum, obtained from the 
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FIG. 5: (Color online) Same as Fig.|4]but for localized initial states: 
Top panel, |4_)|| )||4_)|, two localized monomers; bottom panel, 
|4_)||4+)||4_)|, same plus a localized trimer. 



full BHM. Hole defects with quasi-momenta at the center of 
the band can easily leave the cluster after just a few scatter- 
ing events, Fig. [4] Hole defects prepared at the edges of the 
band remain trapped in the cluster. An additional particle de- 
fect, which itself can not leave the cluster, induces fast quasi- 
momentum redistribution of the hole defects, large fraction of 
which can now leave the cluster. 

The same effect is observed with localized defects. Fig. [5] 
For hole defects alone, about one third of their popula- 
tion leaves the cluster (note that the localized initial state 
of each defect has uniform distribution of quasi-momentum 
k e [— 7r,7r], and not energy Ej^ G [—4 J, 4 J]), while an ad- 
ditional localized particle defect increases this fraction signif- 
icantly. 

In Fig. [6ja) and (b) we show the time evolution of the total 
population of monomers outside the dimer cluster pertaining 
to the cases illustrated in Figs. |4] and [5j respectively. Again, 
hole defects with quasi-momenta in the center of the band eas- 
ily escape the cluster even without the assistance of a particle 
defect, Fig.[6ja). Conversely, for the hole defects with quasi- 
momenta at the edges of the Bloch band in the cluster, very 
little population is found outside the cluster in the long time 
limit (the small fraction of monomer population in the vacuum 
is due to the finite binding energy U of the dimers). Adding 
a particle defect in the cluster significantly increases the frac- 
tion of monomers outside the cluster; we find that the increase 
is always larger for a particle defect in the center of the band 
than for a localized one. 

For the initially localized hole defects, Fig.[6jb), and with- 
out assistance of a particle defect, we find that, as expected, 
about a third of their total population occupying the center 
of the Bloch band leaves the cluster in the long-time limit. 
A particle defect can further increase the portion of escaping 
population of the hole defects by redistributing their quasi- 
momenta over the entire band. 

Note that the results of numerical simulations for the sys- 
tem with a particle defect are reliable for shorter times as com- 
pared to the simulations with the hole defects only, which is 



due to the larger entanglement created dynamically upon the 
trimer-monomer collisions. 

So far we have been restricted to the treatment of only two 
monomers and one trimer and for relatively short times, be- 
cause in the full BHM the fast growing entanglement in the 
system limits the numerical method. With the effective model 
containing only the hole and particles defects, we can simulate 
the dynamics for much longer times with the same numerical 
accuracy, as can be seen in Fig. [6] The perfect agreement be- 
tween the full and effective models allows us to employ the 
effective model for simulating larger systems and for longer 
times. 

Fig. [7] shows numerical results for a system containing ini- 
tially up to four defects. As expected, the evaporation works 
for the larger systems as well. Most importantly, in the pres- 
ence of a particle defects, the number of hole defects left in 
the cluster in the long-time limit falls well below unity (ex- 
trapolating the curves to somewhat larger times than shown in 
Fig.[7] if necessary). 

VI. TWO SPECIES BOSE-HUBBARD MODEL 

We have seen in the previous sections that, in the single 
species BHM, the hopping amplitudes of a monomer inside an 
n = 2 MI cluster and on an empty lattice differ by a fixed fac- 
tor of 2. More flexibility is offered by the two species BHM, 




FIG. 6: (Color online) Total population (integrated density, 
X^jii (^j^j) + X^j=58 (^j^J )) of monomers outside the dimer clus- 
ter, (a) Initial states of the cluster are: in the upper black branch, | — 
V2)M-)iK/2)8"_ (solid line), | -_^/2)M4+)iK/2)^ (dashed 
line), and \-7i/2)l \7i/2)l'^\7v/2)l (dot-dashed line); in the lower 
blue branch, Ml' \-)l\0)l- (solid line), Ivr)^" |4+)i|0)8" (dashed 
line), and |7r)8" |7r/2)8+|0)8" (dot-dashed line), (b) Initial states 
of the cluster are: |4_)i|-)i|4_)i (solid line), |4_)i|4+)i|4_)i 
(dashed line), and |4_)i|7r/2)8+|4_)i (dot-dashed line). The inter- 
action strength is [/ = 100 J. Numerical parameters are the same as 
in Figs.|4]and[5] and the curves terminate when the accumulated cut- 
off error equals 10~^. The gray lines are obtained from the equiva- 
lent effective model, with the time step increased to 1/10 J. 
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FIG. 7: (Color online) Total population + 
5]]^^°gg(a]aj)) of monomers outside the dimer cluster 
of 32 sites surrounded by empty lattice, | )64» on both 
sides. (a) Initial states of the cluster are: in the upper 
black branch, | - 7r/2)^- |.)i|7r/2)^- 1 - 7r/2)^- (solid 
line), I - 7r/2)^-|4+)i|7r/2)^-| - 7r/2)^- (dashed line). 



and 
line); 



- ^/2>M7r/2>|+|7r/2)M - 7r/2>j 
in the lower blue branch. 



(dot-dashed 



k>ri->8io>rk>r 

(solid line), |7r)|" |4+)^|0>8" |7r)|" (dashed line), and 

|7r>M7r/2>^+|0>M7r)~ 



(dot-dashed line). 



(b) Initial 



states of the cluster are: |4_)8|-)8|4_)8|4_)i (solid line), 
|4_)i|4+)i|4_)i|4_)i (dashed line), and |4_)i|^/2)^+|4_)i|4_)i 
(dot-dashed line). (c) Initial states of the cluster are: 
K)8"|-)i|0)8"]4_)i (solid line), K)M4+)i|0)M4_)i (dashed 
line), and I tt) 8 |7r/2)8+|4_)i|4_)i (dot-dashed line). Simulations 
were performed with the effective model. Bond dimensions x — 300 
are used, and the time step size is 1/10 J. The curves terminate 
when the accumulated cut-off error equals 10~^. 



which we now briefly discuss. The Hamiltonian for the sys- 
tem is 

H = - Ja ^(a]a,+i + H.a.) - Jb + H.a.) 

j 3 



2 ' 2 

3 3 



(14) 



where cij (bj) are the bosonic operators for the particles of 
type a (b) hopping between adjacent sites with the rate Ja 



( Jb), while Ua^U\) and Uah are the intra- and inter- species on- 
site interactions. 

Assuming the conditions U^, /7b, U^h, \Ua^U\)- 2Uah\ > 
JajJh, we first consider the situation where each lattice site 
is either empty or contains a single a-b dimer, that is, a pair 
of strongly interacting (via Ug,h) particles a and b localized on 
the same site. Upon adiabatic elimination of the non resonant 
states containing unpaired particles on neighboring sites 1 121 . 
we obtain an effective Hamiltonian of the form of Eq. (|2]), 
where now the dimer hopping and nearest-neighbor interac- 
tion are given by 



J 



2JaJb 



B 



ab 



2Jl 



J' 



ab 



(15) 



With all the interactions repulsive, the anisotropy parameter 

A = s/2J=:^fi + ^V^^fi + ??') (16) 



is larger than 1 for any finite t/a/^^ab or Ui^/Uah, and the MI 
cluster of a-b dimers is stable. But for f/a/^^ab, Uh/U^h 
oo, corresponding to the band insulator for two fermionic 
species, A = 1 and the dimer cluster is unstable. 

Inside the n = na + nb = 2 (rig, = nb = 1) MI clus- 
ter, a hole defect of type a (unpaired particle b) is created by 
dj, see Fig. [ije). The defect hops in the cluster with the rate 
Ja while outside the cluster its hopping rate is Jb. It must 
be stable and not resonantly converted into a pair of particles 
b and a single b-hole (unpaired particle a), which requires 
that U\y — /7ab ^ Ja^Jh- Neglecting the second-order cor- 
rections of the order of Jf b/^a,b,ab, we have the effective 
single-particle Hamiltonian Q with Ja = Ja and Jb = Jb- 
Using the results of Appendix |A| we calculate the transmis- 
sion probability T{k) of the particle through the domain wall 
separating the regions A and B for various Jb/Ja, which is 
shown in Fig. |2] At Ja = Jb (o^ = 1) we find an almost per- 
fect transmission for all /c, up to a small correction due to finite 
interactions. The above results equally apply to a hole defect 
of type b (unpaired particle a) with the replacement a ^ b. 

We have performed numerical simulations of the dynamics 
of several defects in a dimer cluster surrounded by vacuum 



using the full model of Eq. ( 14). For computational reasons, 
we truncate the local Hilbert space to three bosons of each 
species per site, which is justified by the facts that, due to the 
strong interactions, the occupation of a single site by more 
particles can safely be neglected. 

In Fig. [8] we show the behavior of two unpaired particles b, 
or a-holes, moving in the cluster with different initial veloci- 
ties. In the case of Ja = Jb (top panel), both defects almost 
completely leave the cluster upon the first encounter with its 
walls. For Ja ^ Jb, only partial transmission of each defects 
is recorded, which depends on its initial quasi-momentum, as 
per Fig. [2] As an example, at Ja = 2 Jb (central panel) the 
unpaired particle b with /c = 7r/2 can leave the cluster, while 
that with k = tt/A can not, as its quasi-momentum is close to 
the lower band edge. 

Fig. |9] illustrates the results for a pair of initially localized 
defects of different type. Again, for Ja = Jb, both defects 
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FIG. 8: (Color online) Density of unpaired particles b, or a-holes, 
in the cluster (left column), and particles a (right column), in the 
lattice with a MI cluster of a-b dimers spanning 24 sites surrounded 
by empty lattice, | )§2. on both sides. The initial state of the cluster 
I — 7r/2)g~^|-)||7r/4)g~^ corresponds to two a-hole defects moving 
to the left with velocity 2Ja, and to the right with velocity Ja, 
respectively, while all particles a are dimerized with particles b (no 
b-hole defects). The parameters are Ua, — Uh — 60 Jb, Ush — 40 Jb, 
and Ja = Jb (top panel), Ja = 2Jb (central panel), and Ja = \Jh 
(bottom panel). A TEBD algorithm with bond dimension x = 100 is 
used for the time evolution with a fourth order Trotter decomposition 
and time step size 1/50 Jb, with the particle number conservation for 
each species explicitly included in the MPS. 



easily leave the cluster through its walls, but when Ja / Jb, 
only a fraction of the population of each defect leaves the 
cluster after the first collision with its wall. Note, however, 
that since the two types of hole defects have different effec- 
tive mass, their collisions with each other and the walls of the 
cluster can effectively redistribute their quasi-momenta, and 
no trimer defects are required to purify the MI cluster. 



FIG. 9: (Color online) Dynamics of an initially localized unpaired 
particle b, or a-hole, in the cluster (left column), and an unpaired 
particle a, b-hole, in the cluster (right column), for the initial cluster 
state |4_a)8| )8|4-b)8- All parameters are as in Fig.js] and the bond 
dimension of the TEBD is x = 200. 



on top of the n = 2 MI cluster and on an empty lattice. 

In the case of dimers composed of two different bosonic 
species, the defect evaporation proceeds by itself, without the 
need of any catalyzing species. 

The system studied in this paper is amenable to experimen- 
tal investigations with cold atoms in optical lattices |3|. To 
prepare the cluster of dimers surrounded by lattice vacuum, 
one starts with an optical lattice superimposed by a shallow 
confining potential populated by the MI phases with occu- 
pation numbers of n = 0,1,2 in successive spatial shells 
CniEIl, followed by removal of all the atoms outside the cen- 
tral n = 2 MI region 1 22 1 . The homogeneous lattice potential 
is then achieved by turning off the shallow confining potential, 
while the spatial distribution of the defects and their dynam- 
ics inside and outside the dimer cluster can be resolved using 
non-destructive single-site addressing techniques [>23fi25 J . 



VII. SUMMARY 

To conclude, in one-dimensional MI clusters of repulsively- 
bound dimers of bosons | 9 |, hole defects (unpaired particles, 
or monomers) can evaporate through the cluster boundaries, 
taking away the entropy of the system. In the case of dimers 
of identical bosons, only part of the monomer population can 
leave the cluster unassisted. Complete evaporation of the hole 
defects is possible in the presence of catalyzing particle de- 
fects (trimers), which efficiently thermalize the hole defects 
via the quasi-momentum redistributing collisions. The par- 
ticle defects themselves can not leave the cluster, due to the 
large energy mismatch 2/7 between a single excess particle 
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Appendix A: Transmission of a particle through a domain wall 

Here we calculate the probability of transmission of a par- 
ticle with quasi-momentum k through a domain wall, as per 
Eq. ([5]). For the particle incident from the left, we solve the 



9 



stationary Schrodinger equation using the standard scattering 
ansatz for the wave function, 



j>0 



(Al) 



where p and r are the complex reflection and transmis- 
sion ampHtudes. The energy eigenvalue is = e\^^ = 

— 2Ja cos(A:) = — 2Jb cos{k') = Ej^, \ and therefore the re- 
fraction is given by 



k' = cos ^ 



cos(/c) 



(A2) 



Thus the transmission vanishes if cos(/c) > a, where a = 
Jb/ Ja- (In this paper we are primarily concerned with the 
case of a = 1/2, except for Sec.|Vl]) 

Continuity at j = implies 1 + p = r, which together with 
the Schrodinger equation at j = 0, 



yields 



Jee*^' + Jac - 2 Ja cos(/c) 



(A3) 



(A4) 



by 



— Jbg*^' — Jag*^ + 2 Ja cos(/c) 
The current density in the two parts of the system is given 



One can readily verify that 

For the state of Eq. ( |Al| i, we have 

/in = 2JAsin(A;), 
/ref = -2Ja sin(/i;)/j>, 
/trans = 2 Jb sin(fc')T*T, 



/.-l)- 

(A6) 



(A7a) 
(A7b) 
(A7c) 



so that the reflection and transmission probabilities are 

/ref 



R = 
T = 



fi 

^trans 



fir 



- P 

sin(A:') ^ 

= T T. 

sin(/c) 



(A8) 
(A9) 



On can verify that T + = 1 as it should. 



Appendix B: Two particle Hamiltonian in quasi-momentum 
space 

We consider a pair of distinguishable, locally interacting 
particles on a lattice described by Hamiltonian 



H 




H.a. 



(Bl) 



where the periodic and open boundary conditions correspond, 
respectively, to 7 = 1 and 7 = 0. The operators and cre- 
ate soft-core particles interacting via U, which is convenient 



for the exact numerical simulations presented in Sec. |IVB 

In quasi-momentum representation, k = 27ny/L {y = 
L-f + IJ . . . Lf J), we have al = ^ Ylf^^ e'^^a] and 
£^ = Xlj^i ^iid two such particles have a proba- 

bility to be on the same real lattice site. The Hamiltonian 
then reads 



H 



-2Ja y^cos(A:)a^afc 

k 

-2Jt ^cos(A:)44 4 



(l-7)^a 

L 



) 



L 



Y 2^ CLk^k'ik"'^{k+k" 



(B2) 



Appendix C: Effective theory for monomers and trimers 

The non-locality of the effective theory presented in 
Sec. |VA| might seem surprising at first sight. From the point 
of view of quantum information theory, however, the Hamil- 
tonian ([9]) is still local, in the sense that the commutator 
[Hj^Hjf] vanishes except for / = j ± 1, despite the fact 
that the support of any two Hj^Hj' has a large overlap. This 
property should always be conserved in any effective theory, 
since it guarantees that correlations in the model system travel 
with the same maximal velocity as in the full system |26, 27 1. 
This property also permits the application of the TEBD nu- 
merical method, in conjunction with the conservation of the 
total particle number, to the effective model. For then all the 
basis states used in the TEBD (eigenstates of the reduced den- 
sity matrices for all bi-partitions of the lattice) are, by con- 
struction, the eigenstates of the total particle number in the 
corresponding subsystem. Since the total particle number is 
the only observable that enters Hamiltonian d9l via P^^ , this 
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FIG. 10: (Color online) Comparison of the full Bose-Hubbard dy- 
namics, Eg. (Tl with U = lOOJ, (a) and (b), with the effective 
model, Eq. (|9]), (c). The initial state contains a MI cluster of dimers 
on sites j > I and localized monomers at sites j = — 1 and j = 2. 
The density of dimers is shown in (a), and the density of monomers 
in (b)-(d). Note that the cluster boundary is shifted upon parti- 
cle crossing, which manifests in (a) as a smoothing of step in the 
dimer density. The effective Hamiltonian without moving bound- 
aries, H ^Y.^. 



-1 M'^m 

:1 



yields the dynamics of (d). 



type of non-locality does not introduce additional difficulties 
in the use of the TEBD method. 

The effective model can also be extended to higher orders 
in perturbation theory. In second order, this introduces nearest 
neighbor interactions, local potentials and effective exchange 
between monomers and trimers. All these terms are of the or- 
der of J'^ /U and depend on 6, which can now assume four 
different values depending on the type of bond between sites 
j and j + 1. Another term of the same order describes hole 
defect hopping to the next-nearest neighbor site in the clus- 
ter. As this is spanning three sites, it also depends on the state 
of the central site and requites more values of 9. The pres- 
ence of such a longer-range term would necessitate a more 
general numerical simulation algorithm than TEBD. We have 
verified, however, that the effective Hamiltonian ([TO]) contain- 
ing only the terms first order in J already captures all the es- 
sential physics discussed in this paper. 

In Eig. [T0| we compare the dynamics of hole defects ob- 



tained from the full and effective models, which agree very 
well for large interaction strength U ^ J. Observe, how- 
ever, that a local theory neglecting the motion of the cluster 
boundaries, Eig. [TOjd), and therefore violating the conserva- 
tion of the total number of dimers and bare particles, does not 
describe the dynamics quantitatively correctly. 



Appendix D: MPS representation of the initial state 

Here we show how to construct an exact MPS representa- 
tion for a lattice containing fixed number of bosons each in a 
certain single particle eigenstate. The resulting MPS will be in 
the canonical representation |28| and symmetric |29|, i.e., it 
will be an eigenstate of the total particle number by construc- 
tion. The construction is analogous to that of matrix product 



operators for fixed total particle number II301 . 

The single particle state is given by a normalized wave- 
function In the examples of Sec. VB we have = 
_^g*fej with fixed quasi-momentum k. The corresponding 

state 1^^) = ZljLi can be written as 

1^') = (V^^A + - ^m4) |0)a ^ |0)b, (D1) 

where sub-lattice A spans sites I to m and sub-lattice B is 
from m + 1 to L, while qm = Xljli ^j^j the single particle 

LllLJ' U'V^lllg, 111 J~V. X IIV^ UWOWlllV^ V^lV^CltlWll WpV^l CltWl O U/^ 

and ttg are defined by 



probability of being in A. The bosonic creation operators 

L 



a 



I 



a/I - Qrr 



j=m+l 



(D2) 

The state of the lattice with N particles in the same single 
particle state can then be expressed as 



1^ 



N\ 



;^/=(v^^A + \/l - qma'l)^\0)A <^ |0)b 
x(\/9^aA)'(\/l - <lma\i)^~^\0)A <S) |0)b, 



and the density matrix of the system is 



1 ^ 
~ AH E 



l,l'=0 



x(v^^a) |0)a(0|a(v^^a) 
^{V^ - gm^e)^ ^|0)b(0|b(\/1 - QmaB] 
The density matrix of subsystem A is 
PA = TrB[|*^)(*^|] 



N-l' 



(D4) 



I ^ 



N 



1=0 ^ 



{N-l)\{l-qm) 



N-l 



X (v^^a)^|0)a(0|a(v^^a)'. 



(D5) 



Note that pA has at most x = ^ + 1 nonzero eigenvalues, 
one for each possible distribution of the N particles between 
A and B. With P^^^ the projector onto the / particle sector of 
subsystem A, the probability of finding / particles in A is 

TrA[pAP/^l] = ^f:(^)'(iV-0!^!(l-9n^)^-'(g™)' 

N 



/=0 



(D6) 
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which is a binomial distribution. 

We can now construct |^^) as a matrix product state in 
the canonical ||28]| form. Given a bi-partition of the lattice, its 
Schmidt decomposition is 



N 



1^ 



N 



1^ 



N-l\ 



(D7) 



1=0 



with \¥)a = ^{aiy\0)A and |^^-% = 

/ ^ (^r)^^~^^|Q)b. The MPS will have bond di- 

mension of x = ^ + 1- The probability of finding / particles 
to the left of bond m is 



(D8) 



We then continue with the Schmidt decomposition at the fol- 
lowing bond. The remaining task is to determine the coeffi- 
cients of 



N N 



1=0 r=l 

x|*')a®|*'-')„^+i®|*^-'")b'. (D9) 
The A tensors are already known from (|D8|). The sub-chain 



comprises sites m + 2 to L. Thus (A 



■ Ir 



1+1] 1 2 



is the probability of finding N—r particles on the right of bond 
m + 1 and I particles on the left of bond m, resulting in 



|-p[m+l] |2 



B_ 



-{l\r) 



r\{N-l)\ 
{r-l)\N\ 



^r-l 



m+1 



(1 



(DIO) 



J-N 



For the phase to be correct, we obviously have to set 



[m+l]N 



(r - /) arg (0m+i) 



(DID 



Equations ( |D8| ), ( DIP ), and ( |D11| ) completely determine the 
tensors T and A. Note that in this particular case, the value of 
the bond index of A^^^ has a physical meaning of the number 
of particles to the left of bond m. The resulting MPS is an 
eigenstate of the total particle number, which can be used in 
TEBD implementations that take advantage of particle num- 
ber conservation explicitly, as in this paper. 



The construction is more complicated when one intends to 
prepare particles in different single particle states a = 
1, 2, . . . , M. From simple combinatorial considerations, we 
deduce that the Schmidt rank will be x = Ylai^ + ^a)^ i-^., 
exponentially large in the number M of different single parti- 
cle states. (This implies that, as a starting point for dynamical 
simulations, one can construct an exact MPS for the ground 
state of non-interacting bosons, as done in [31 1, but not for 
non-interacting fermions.) The exact expression in terms of 
the qm,a will contain overlaps between the different single 
particle states, which in general are finite in any subsystem 
even if the single particle states are orthogonal on the entire 
lattice. 
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